Probing the parameter space of HD 49933: a 
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Abstract. We present two independent methods for studying the global stellar parameter 
space (mass M, age, chemical composition Xo, Zq) of HD 49933 with seismic data. Using a 
local minimization and an MCMC algorithm, we obtain consistent results for the determination 
of the stellar properties: M = 1.1-1.2 M©, Age ~ 3.0 Gyr, Z ~ 0.008. A description of 
the error ellipses can be defined using Singular Value Decomposition techniques, and this is 
validated by comparing the errors with those from the MCMC method. 



1. Introduction 

HD 49933 is a main sequence solar-type star that was observed by CoRoT. It is the first stellar 
object where solar-like oscillations were clearly detected in the photometric signal [HE]. There 
has been some controversy over the original labelling of the oscillation frequencies with their 
mode degrees, and in this work we present results based on the data published by [lj, although 
we note that since this publication there has been a preference among the scientific community 
towards an inverted mode-labelling [2J. We present two independent methods of fitting the 
observational data, while placing an emphasis on defining the boundaries of the parameter space 
where the model of this star lies. The best-fitting models are determined using the Levenberg- 
Marquardt (LM) and Markov Chain Monte Carlo (MCMC) algorithms, while the uncertainties 
and the form of the parameter space that is constrained by the set of observations is described 
by both Singular Value Decomposition (SVD) and MCMC. 

We use the frequencies from [1], and we compute the frequency separations to use as the 
observational data to model. The stellar models are computed using ASTEC [3] which include 
the EFF equation of state, OPAL 95 opacities together with Kurucz low-T opacities, NACRE 
reaction rates and overshooting, and the oscillation frequencies are calculated with ADIPLS [4J. 

2. Local minimization of stellar observables 

Using as observational constraints the £ = frequency separations, T e ff,Lo,[M/H] from [lj, 
we obtain the following set of stellar parameters for HD 49933 using LM: M = 1.12 M , Age 
= 2.9 Gyr, Z = 0.007, a = 1.80, and R = 1.39 R . X = 0.70 was fixed, and the analysis 
was repeated using various initial parameter values [5j. In Fig. [I] each dot is the result of a 
minimization, clearly showing a sensitivity to the initial guess. This sensitivity is not uncommon 



for local methods, however, when we repeated our analysis using the data from [2] we found a 
rather stable solution of M = 1.20 M , Age = 2.5 Gyr, Z = 0.0105, a = 1.70 and R = 1.44 R . 

3. Global stellar parameter solution 

Markov Chain Monte Carlo (MCMC) algorithms allows one to perform stochastic samplings of 
probability densities using properties of Markov Chains. This is a Bayesian methodology and 
we use it to estimate the posterior Probability Density Function (PDF) of the parameter (s) 
of our stellar model. The strength of this approach is that, once the PDF is derived, we can 
apply the classical tools of statistical inference to estimate the stellar parameters and associated 
confidence intervals — a long-standing problem in stellar physics [6j [71 |5j. This methodology 
was applied successfully to a Cen A [8] . 

We ran the MCMC algorithm using two different values of initial parameters (M = 1.1, 1.2 
Mq, Zq = 0.007,0.010) while including a prior on Zq for the first run only. Using the same 
observational constraints as those described in Sect. 2, we computed the marginal distributions 
of the global stellar properties. Using their mean values as the best-fitting parameter and 
considering the standard deviation of the sample as the corresponding 1-a error bar, we obtain 
M = 1.104 ± 0.010 M , Age = 3.16 ± O.llGyr, Z = 0.0065 ± 0.0006, and a = 2.002 ± 0.080. 
However, we also find that the second MCMC run hints at a possible (much less populated) 
alternative solution with a mass near 1.2 M . 

4. SVD describing the error ellipses 

When we obtain the best-fitting set of parameters P, then locally the models can be described 
as linear. With this assumption, we proceed to calculate the SVD (D = UWV T ) of the matrix 
D — f9 a_1 ' wnere O are the observables and a the observational errors. The matrix V describes 
the TV-dimensional (in our case N = 4) parameter correlations (essentially 4 4-element vectors), 
while W describes the magnitudes of these vectors. Fig. [2] shows the (longest) two-dimensional 
projections of the four-dimensional vectors, for the parameters of M, age and Zq. These 2-D 
vectors are represented by the red arrows in the figures, and the black ellipses are the error- 
ellipses defined by these vectors. 

5. Uncertainties 

Fig. [2] shows the error ellipses described by SVD. The models generated from the MCMC 
runs are represented by the blue dots, with the size of the dots proportional to the likelihood 
values calculated from the MCMC (i.e. larger dots are more likely to be the true models). 
The error ellipses clearly encloses most of the (good) models from the MCMC method. This 
is strong evidence that, although we are assuming linearity for these models, the analytical 
approximation describes quite well the correlations that we expect to find, and the magnitudes 
of these correlations. 

Because SVD gives a good indication of the shape of the parameter spaces, we can use it to 
investigate the effects of including extra observational constraints. In Fig. [3] we compare the 
two-dimensional error-ellipses for various parameters when using the same constraints as Fig. [2] 
(black ellipses), and then including the large frequency separations from the t — 1 modes (blue 
ellipses). The allowed parameter spaces defined by the extra set of constraints is reduced in 
volume by a factor of 2 for the 4 parameters, and the individual uncertainties are reduced by a 
significant amount for Z§ and a. 
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Figure 1. Levenberg-Marquardt results for the best-fitting parameters using the £ — large 
frequency separations and T e ff , L , [M/H] data from [lj as constraints on the models. The 
various colours indicate different initial guesses of the parameters, all having fixed Xq = 0.70 
except for the red dots which have fixed Xq = 0.73. Each colour has a set of solutions, because 
we also varied the initial guess of the mass. 
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Figure 2. 2-D probability distributions for M and Age (left) and M and Zo (right) from the 
MCMC run. The size of the blue dots is proportional to the likelihood of the corresponding 
model. The red arrows are the two-dimensional scaled projections of M and Age as defined by 
SVD, while the black ellipse is the analytical error-ellipse given by these vectors. 
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Figure 3. Error ellipses for M and Age (left), M and Zq (centre) and Age and a (right) 
as defined by SVD, assuming a set of constraints comprising T e ff,L and the large frequency 
separations for £ = only (black), and then including the £ = 1 frequency separations (blue). 



6. Conclusions 

We have analysed the seismic data from HD 49933 and we have found that the global stellar 
properties using two independent methods (local and global) are in agreement. Using the 
mode-identification and frequencies given by p] and using the frequency separations as the 
observational data to compare to the models we obtain M — 1.12 M ± 2%, Age = 2.9 Gyr 
± 5%, X = 0.70 (fixed), Z = 0.007 ± 18%, a = 1.80 ± 9%, and R = 1.39 R from the 
local minimization method. This solution is supported by the results from the global method: 
M = 1.1 M , Age = 3.16 Gyr, Z = 0.0065 and a = 2.00. The MCMC run also finds a possible 
alternative solution with M = 1.2 Mq. A similar solution is also found from the local method 
while using the reversed mode-tagging [2]: M = 1.20 M , Age = 2.5 Gyr, Z = 0.0105, a = 1.70 
and R= 1.44 R . 

We also studied the parameter space defining the correlations and uncertainties using the 
analytical approach of SVD. The results from the global MCMC validated the analytical 
description. In particular we found that including the £ = 1 frequency separations should 
reduce the volume of parameter space by at least a factor of two in all dimensions, with the 
individual parameter uncertainties reducing significantly for Zq and a. 
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